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Abstract 



We propose a self-consistent method for electronic structure calcula- 
tions of correlated systems that combines the local density approximation 
(LSDA) and the dynamical mean field theory (DMFT). The LSDA part is 
based on the exact muffin-tin orbitals (EMTO) approach, meanwhile the 
DMFT uses a perturbation scheme that includes the T-matrix with fluc- 
tuation exchange (FLEX) approximation. The current LSDA+DMFT 
implementation fulfills both self-energy and charge self-consistency re- 
quirements. We present results on the electronic structure calculations 
for bulk 3d transition metals (Cr, Fe and Ni) and for Fe/Cr magnetic 
multilayers. The latter demonstrates the importance of the correlation 
effects for the properties of magnetic heterostructures. 
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I. INTRODUCTION 



In the ab initio description of the electronic properties of materials the most widely 
used methods are based on the density functional theory (DFT) [1] implemented within 
the local spin density approximation (LSDA) [2,3] to the exchange and correlation en- 
ergy. Ground state properties of the most of metals, semiconductors, ionic compounds, 
etc., are quantitatively well described by the DFT-LSDA approach. Attempts to ap- 
ply these first principles methods to strongly correlated systems, however, encountered 
many fundamental difficulties [4-6]. Even for elemental transition metals, such as Mn, 
Fe, or Ni, the impact of the correlation effects on the electronic structure turns out 
to be essential [7]. Therefore, one of the most challenging problems in the physics of 
transition metals, their alloys and compounds is to develop simple and efficient elec- 
tronic structure methods that go beyond the LSDA by including important many-body 
effects. 

It has proved a fruitful approach to combine the simple Hubbard model with the 
LSDA technique, providing a DFT scheme "beyond LSDA" [4-7]. Unfortunately, the 
simplest realization of such an approach, the LSDA+U scheme [4], cannot describe the 
many-body effects beyond the Hartree-Fock approximation. These effects are connected 
with the frequency dependence of the electron self-energy. In order to include dynam- 
ical effects the LSDA+U scheme was combined with the dynamical mean field theory 
(DMFT) [5,6]. The DMFT maps lattice models onto quantum impurity models subject 
to a self-consistent condition in such a way that many-body problem for crystal splits 
into one-body impurity problem for the crystal and many-body problem for an effective 
atom. In fact, the DMFT, due to numerical and analytical techniques developed to 
solve the effective impurity problem [8] , is a very efficient and extensively used approx- 
imation for energy dependent self-energy H(ou). The emerged LSDA+DMFT method 
can be used for calculating a large number of systems with different strength of the 
electronic correlations [7,9,10]. To underline the importance of complete LSDA+DMFT 
self-consistency we mention that the first successful attempt to combine the DMFT with 
LSDA charge self-consistency gave an important insight into a long-standing problem 
of phase diagram and localization in f-electron systems [11]. 

To incorporate the dynamical mean field approach in the band structure calcula- 
tion we adopt the exact muffin-tin orbitals (EMTO) density functional method. The 
EMTO theory can be considered as a screened Korringa-Kohn-Rostoker (KKR) muffin- 
tin method, where large overlapping potential spheres are used for accurate represen- 
tation of the LSDA one-electron potential. A comprehensive description of the EMTO 
theory and its implementation within the LSDA may be found in Refs. [12] and [13,14], 
respectively. 

The paper is organized as follows: Section II presents a general formulation of the 
combined multiple scattering and dynamical mean field approach. The calculation 
scheme from Section III illustrates the multiple scattering solution of the LSDA prob- 
lem via the EMTO method and the many-body solution of the DMFT problem via 
the T-matrix FLEX approach. First-principles results obtained from EMTO-DMFT 
calculations are discussed in Section IV. The paper is summarized in Section V. 
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II. FORMULATION OF THE PROBLEM 



The density functional theory reformulates the N electron problem into one electron 
problem by considering a non-interacting system, where each electron feels an effective 
potential ^e//( r ) created by the rest of the electrons and external fields. Thus, within 
the DFT the solution of the original inhomogeneous system is constructed from the 
one-electron Kohn-Sham equations [2] 

■V 2 + ^(r)] ¥*(e,r) = e^(e,r), (1) 

where a stands for spin. The many-body part of the effective potential /i£ c (r) is an 
unknown functional of the spin densities n a (r) = J2 e \^ a { e > r )| 2 - The most commonly 
adopted approach for //£ c (r) is the local spin density approximation (LSDA), where the 
effect of interactions between electrons is taken into account by substituting locally the 
real system by the uniform electron gas with the density equal to the actual density 
at point r. In this paper we will not distinguish between different specific forms of the 
LSDA. 

In order to include the many-body correlation effects beyond the LSDA we substitute 
the Kohn-Sham equation (1) by the quasiparticle equation 



-V 2 + ^(r)] S"( e ,r) + 



+ ]T \Rlma)^ mmtRlml (e)(Rlm'a\^(e)) = e^(e,v), (2) 

Rl mm' 

where R, I and m denote the lattice sites, the orbital and the magnetic quantum num- 
bers, respectively. \Rlma) are localized orthonormal basis functions, e.g. the partial 
waves for the correlated /—channels. In Eq. (2) the correlation effects are treated at the 
DMFT level, where the essential many-body self-energy, S^, /mram , (e), is a local, energy 
dependent and multi-orbital potential. 

Note that Eqs. (1) and (2) are formulated in terms of wave functions. Consequently, 
the DMFT method has already been implemented in several techniques based on wave 
function formalism, like the linear muffin-tin orbital method [5,6,11]. At the same time, 
accurate self-consistent methods for solving the local Kohn-Sham equation (1) in terms 
of Green's function have been developed within the multiple scattering theory [15-18]. 
The main aim of the present work is to include the many-body correlation effects, 
approximated by means of the DMFT, in the above mentioned multiple scattering 
approach. We note that for a general non-local energy dependent potential the multiple 
scattering theory offers a solution known as the optical potential [19]. However, the 
optical potential is far to complicated to be used in realistic computation. Nevertheless, 
it is proved that the non-local potential could be transformed into an one particle 
energy-dependent operator such that it satisfies a similar one-particle equation with 
local and energy independent potential. 
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III. THE CALCULATION SCHEME 



A. The one-electron Green's function 



Within the multiple-scattering formalism, the one-electron Green's function is de- 
fined for an arbitrary complex energy z as 

z + Vl- v\ ff {v)\ G°(v, r', z) = 6(t - r'). (3) 

For most of the applications, e.g. standard KKR or LMTO methods, the LSDA effective 
potential from Eq. (3) is approximated by spherical muffin-tin (MT) wells centered at 
lattice sites R. Within a particular basis set, the one-electron Green's function is 
expressed in terms of the so-called scattering path operator, g R L, S R >L'(z), as wen as the 
regular, Z RL (z,r R ), and irregular, J RL (z,r R ), solutions to the single site scattering 
problem for the cell potential at lattice site R, viz. 



G a,LSBA (rR + R Yr/ + R , z) = J2 Zrl (z, V R )g^,Uz)Z% LI {z, T R ,) - 

L,L> 

-6 rr ^Jrl(^r)Z rl (z,v r ), (4) 

L 

where L = (7, m) with I < l max (usually l max = 3) and i\r = r — R denotes a point 
around site R. The real space representation for the scattering path operator for the 
muffin-tin potential is given by 

9rl S r>l>(z) = [$RR>$LL't a R 1 L (z) - B RLtR , L ,(z)Y l , (5) 

where t RL {z) stands for the single scattering t-matrix and B R l,r' u {z) are the elements 
of the so-called structure constant matrix. 

Unfortunately, the MT based KKR or LMTO methods have limited accuracy. The 
former method uses non-overlapping spherical muffin-tin potentials and constant po- 
tential in the interstitial, while the latter method approximates the system with over- 
lapping atomic sphere and neglects completely the interstitial and the overlap between 
individual spheres. Recent progress in the field of muffin-tin orbitals theory [12] shows 
that the best possible representation of the full potential in terms of spherical wells may 
be obtained by using large overlapping muffin-tin wells with exactly treated overlaps. 
Within this so-called exact muffin-tin orbitals method [12], the scattering path operator 
is calculated as the inverse of the kink matrix defined by 

k rl,r'l'( z ) = ^RR'hL'D RL (z) - Sbl,wv(z), (6) 

where D RL (z) denotes the EMTO logarithmic derivative function [13,14], and 
Srl,r'L'(z) is the slope matrix [12]. 

Since the energy derivative of the kink matrix, K RL R i L i(z), gives the overlap matrix 
for the EMTO basis set [12], the matrix elements of the properly normalized LSDA 
Green's function become [13,14] 



g rl S Sl'( z ) = H 9 r Y S r^l''( z )K RIILII)RIU {z) - 5 RR ,5 LU r RL (z)i (7) 

R"L" 
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where Irl( z ) accounts for the unphysical poles of K RL R /l>(z). In the case of translation 
invariance Eqs. (6) and (7) can be transformed to the reciprocal space, so that the lattice 
index R runs over the atoms in the primitive cell only, and the slope matrix, the kink 
matrix, and the path operator depend on the Bloch wave vector k. In this case the 
total number of states at the Fermi level E F is obtained as 

N(E F ) = ± £ i [ G R f% A L ,(k,z)dkdz, (8) 

Z7F< RL,R'L' J JBZ 

where the energy integral includes the Fermi distribution. The k-integral is performed 
over the first Brillouin zone, while the z-integral is carried out on a complex contour 
that cuts the real axis below the bottom of the valence band and at E F . 



B. DMFT Green's function and effective medium Green's function 



To incorporate the many-body effects into the Green's function technique we start 
with the LSDA Green's function matrix (7) expressed on the EMTO basis set. The 
LSDA+DMFT Green's function, G RL R , L ,(k, z), defined for Bloch vector k and energy 
z, is connected to the one-electron LSDA Green's function trough the Dyson equation 



Grl.r'l'O*-, z) 



^-,ct,LSDA / 1 \ 
" RL,R'L'\^ z ) 



n -1 



^RR'^RL,RL 



,(z). 



(9) 



The local self-energy £#lrz/ (z) in Eq. (9) depends on the so-called effective medium 
or bath Green's function Q RL R > L r(z). This, in turn, is calculated from the k— integrated 
LSDA+DMFT Green's function, G RL R , L ,(z) = j BZ G RL R , Lf (k, z)dk, as 



^RL,R'L'( Z ) 



Grl,R>L'{ Z ) 



-i 



+ 3rr'^rl,rl>( z )- 



(10) 



To find the local self-energy we use a spin polarized T-matrix plus fluctuation exchange 
(SPTF) approximation [20]. The many body problem is solved on the Matsubara 
contour, defined by the fermionic frequencies u n = (2n + l)nT, where n = 0, ±1, 
and T is the temperature. A Pade analytical continuation [21] is used to map between 
the complex energies z, used in the EMTO iterations, and the complex energies iu n , 
corresponding to the Matsubara frequencies and expressed relative to the Fermi level 
E F . Next we describe the solution of the effective impurity problem. 



C. The solution of effective impurity problem 

The many-body problem is solved using the SPTF method proposed in Ref. [20], 
which is a development of the earlier approach [22]. The SPTF approximation is 
a multiband spin-polarized generalization of the fluctuation exchange approximation 
(FLEX) of Bickers and Scalapino, but with a different treatment of particle-hole (PH) 
and particle-particle (PP) channels. Particle-particle (PP) channel is described by a 
T-matrix approach [23] giving a renormalization of the effective interaction. This ef- 
fective interaction is used explicitly in the particle-hole channel. Justifications, further 
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developments and details of this scheme can be found in Ref. [20]. Here we present the 
final expressions for the electron self-energy. The sum over the ladder graphs leads to 
the replacement of the bare electron-electron interaction by the T-matrix which obeys 
the integral equation 

< 13\T™' (in)\24 >=< 13M24 > -T]T ]T G&(iu))G&(iQ. - iu>) < Q8\T aa ' (ity |24 >, 

u> 5678 

(11) 

where the matrix elements of the screened Coulomb interaction, < 13|i>|24 >, are 
expresed using the average Coulomb and exchange energies U, J [20] . In this section, 
for sake of simplicity, we use the short notation 1 = Rim. In the following we write 
the perturbation expansion for the interaction (11). The two contributions to the self- 
energy are obtained by replacing of the bare interaction by a T-matrix in the Hartree 
and Fock terms 

^ n (iuj) = T]T ]T < 13\T aa ' {iQ)\2A > g^{itt - iu) 

tt 34cr' 

E^ F (iuj) = -Tj2J2< 14|T CTCT '(^)|32 > g£{iQ-iuj). (12) 

fl 34ct' 

The four matrix elements of the bare longitudinal susceptibility represents the density- 
density (dd), density- magnetic (dm ), magnetic-density (m°d) and magnetic- magnetic 
channels (m°m°). The matrix elements couples longitudinal magnetic fluctuation with 
density magnetic fluctuation. In this case the particle-hole contribution to the self- 
energy is written in the Fourier transform form 

^(r) = T,W^ 2 (r)g^(r), (13) 

34<t' 

r being the imaginary time. The particle-hole fluctuation potential matrix W aa ' (icu) is 
defined in FLEX approximation [24,22] 

We emphasize that all of the above expressions for the self-energy, in the spirit of the 
DMFT approach, involve the Weiss (bath) Green function (10). The total self-energy 
is obtained from Eqs. (12) and (13) 

^(iuj) = £ CT ' TH M + £ CT ' TF M + £ CT ' PH M. (15) 

Since the LSDA Green's function already contains the average electron-electron inter- 
action, in Eqs. (9) and (10) the static part of the self-energy S CT (0), is not included, i.e. 
we have 

t a {iu) = Y?{ioj) - S a (0). (16) 
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D. The charge self-consistency loop within the EMTO-DMFT scheme 

After the self-energy ^rl,rl'( z ) * s determined as the self-consistent solution of 
the effective impurity problem, the many-body LSDA+DMFT Green's function, 
G a RL R , L ,(k,z), is calculated using Eq. (9). The EMTO-DMFT number of states at the 
Fermi level is given by the multi center expression (8), written for the LSDA+DMFT 
Green's function. 

The charge and spin densities in the EMTO formalism are represented in one center 
form around each lattice site R, i.e. 

n°(v)=^n RL (r R )Y L (r R ), (17) 

RL 

where Yl(?"r) are the real harmonics. Inside the Wigner-Seitz cell the partial compo- 
nents n RL (r R ) are expressed in terms of the density matrix T> RL , L (z) as 

4W = ^T, I E C L L"L<Z a m »{z, r R ) V RL „ L ,(z) Z RV (z, r R )dz, (18) 

ZTll JEp j^ii^i 

where C^„ L , are the real harmonic Gaunt coefficients. The density matrix is obtained 
from the path operator as described in, e.g., Ref. [14]. Within the present EMTO- 
DMFT scheme the LSDA+DMFT path operator, g RL ^/(k, z), is determined according 
to Eq. (7) using the LSDA+DMFT Green's function and the LSDA overlap matrix. 
Here we implicitly make the assumption that the LSDA+DMFT Green's function can 
be expanded on the same basis set as the one-electron Green's function. In other words, 
instead of the solutions of Eq. (2) for a single scatterer, we use the LSDA single-site 
solutions of Eq. (1) for a single scatterer to express the LSDA+DMFT Green's function 
in the same form, as we express the one-electron Green's function in Eq. (4). 

Finally, for charge self-consistent calculation we construct the new LSDA effective 
potential from the spin and charge densities, n(r) = n T (r) + n^(r). The Poisson's 
equation is solved using the spherical cell approximation [14], and the exchange and 
correlation term is calculated within the LSDA. 

IV. RESULTS AND DISCUSSION 

The role of the correlation effects in the electronic structure of 3d transition metal 
series is far from being completely understood [7]. In general, spin-polarized band 
structure calculations give adequate description of the ferromagnetic ground state for 
the most of metals. At the same time, there are obvious evidences of essentially many- 
body features in photoemission spectra of Fe [22], Co [25], and Ni [26]. Few examples are 
the 6 eV satellite in Ni density of states, broadening of the angle-resolved photoelectron 
spectroscopy (ARPES) features due to quasiparticle damping, narrowing of the d-band, 
essential change of spin polarization near the Fermi level, etc. Although, there are no 
direct experimental information yet, one can assume that the many-body effects can also 
be important in the case of magnetic multilayers and other heterostructures containing 
transition metals. The importance of correlation effects on transition metal surfaces 
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has already been demonstrated by the STM observation of an orbital Kondo resonance 
in Cr [27]. Here we will present results obtained using the EMTO-DMFT method in 
the case of bulk Ni, Fe and Cr 3d transition metals and for Fe/Cr multilayer structure. 

A. Numerical details 

The self-consistent EMTO-DMFT calculations were carried out for the experimental 
ground state crystal structures, i.e. fee for Ni, and bee for Fe and Cr. The lattice 
parameters were fixed at the experimental values. The studied Fe/Cr multilayer system 
has tetragonal (001) structure with one type of Fe atoms and two different types of Cr 
atoms. This structural setup describes the situation of one Fe layer embedded in few Cr 
layers. The atoms in the tetragonal unit cell were fixed in the ideal positions, and for the 
lattice parameters we used the bulk Cr lattice constant. The LSDA Green's function 
was calculated for 16 complex energy points distributed exponentially on a semi-circular 
contour. The /c-point sampling was performed on a uniform grid in the Brillouin zone. 
For the LSDA energy functional we used the Perdew-Wang parametrization [28] of 
the results by Ceperley and Alder [29]. The DMFT parameters, average Coulomb 
interaction U, exchange energy J and temperature T, used in the present calculation 
are listed in the last three columns of Table I. 

From the self-consistent density of states (DOS) we have determined the magnetic 
moment fi and the electronic specific heat coefficient 7. The latter is given by relation 

1 = 7r 2 k 2 B N(E F )(l + X)/3, (19) 

where N(E F ) is the electronic DOS at the Fermi level, (1 + A) is the mass enhancement 
factor caused by the electron-phonon interaction. This factor in the case of Ni was 
estimated to be 1.24 [30]. The present theoretical results for the self-consistent magnetic 
moments and electronic specific heats, along with the available experimental data are 
listed in Table I. 

B. Ni 

It has been shown that the main peculiarities of the experimental Ni photoemis- 
sion spectra can be understood within the framework of the LSDA+DMFT approach 
[26]. An exact quantum Monte-Carlo (QMC) solution of the effective impurity problem 
[26] gives impressive quantitative agreement between the experimental and computa- 
tional data, both for photoemission spectra and for temperature dependent magnetic 
properties. Here we will show that the perturbative SPTF approach, employed in the 
present EMTO-DMFT method, also reproduces the main correlations effects beyond 
LSDA in Ni, e.g. the narrowing of the band, reduction of the exchange splitting and 
the appearance of the 6 eV satellite. 

The EMTO-DMFT density of states for fee Ni is shown on Fig. 1. For the present 
choise of the average Coulomb interaction U = 3eV, the position of the 6 eV satellite 
is shifted to the lower energy. This shift and the large broadening of the resonance is 
due to the perturbative approach of the solver of the effective impurity problem [20]. 
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According to the present result, at T = 500 K, i.e. at T/T c 0.8, the reduction of 
the exchange splitting relative to the LSDA value is 43%. This is in good agreement 
with 40% estimated from experimental data corresponding to the same temperature. 
Our results are also in agreement with ARPES measurements [31,32]. The EMTO- 
DMFT magnetic moment of 0.42/xs, see Table I, represents a reduction of 30% from 
the LSDA value. This reduction is comparable with 25% evidenced from experimental 
magnetic moments. 

Apart from the many-body self-consistency, the present implementation fulfills also 
the charge self-consistency. This allows us to monitor the effect of DMFT on the 
LSDA charge and magnetic moment densities in the real space. Fig. 2 shows the LSDA 
and LSDA+DMFT magnetic moment densities in Ni along the (110) direction, which 
corresponds to the nearest-neighbor distance in the fee unit cell. Comparing the two 
densities one can interpret the reduction of the magnetic moment in the LSDA+DMFT 
approach as a slight narrowing of the real space extension of the d wave functions. 

The energy dependence of the self-energy for Ni, plotted on Fig. 3, near the 
Fermi level shows the typical Fermi liquid behavior. For the imaginary part we have 
—Im oc E 2 , meanwhile the real part of the self-energy has a negative slope 

dRe Y>(E)/dE < 0, where E is the electron energy relative to the Fermi level. 

Within the LSDA for the electronic specific heat coefficient we obtain 5.43 
mJ/K 2 mol, which underestimates the experimental value from Ref. [33] by more than 
20%. This LSDA value is in good agreement with previous calculations [34,35]. On the 
other hand, within the LSDA+DMFT approach the electronic specific heat coefficient 
is 6.78 mJ/K 2 mol, which reduces the discrepancy between theory and experiment by 
19%. 

C. Fe 

In Fig 4 we compare the LSDA and the LSDA+DMFT density of states for bee 
Fe. The present LSDA magnetic moment 2.25/xb is reduced to 2.23/zs within the 
LSDA+DMFT approach. We note that our LSDA magnetic moment is in excellent 
agreement with the one obtained in former ab initio calculation [36] . The electronic spe- 
cific heat coefficient increase from the LSDA value of 2.43 mJ/K 2 mol to 2.53 mJ/K 2 mol 
in the LSDA+DMFT calculation. These values can be compared with the experimental 
data in the range of 3.11 — 3.69 mJ/K 2 mol [37]. We can see that in bee Fe the correlation 
effects are much less pronounced than in fee Ni. This is due to the Fe large spin split- 
ting and the free-structural dip in the density of states [26]. The energy dependence of 
the self-energy of Fe, see Fig. 5, corresponds again to the Fermi liquid behavior, similar 
to the case of Ni. 

The d band exchange splitting in Fe in the LSDA+DMFT is slightly decreased 
in comparison with the LSDA result. The temperature dependence of the exchange 
splitting in Fe was determined by spin-resolved photoemission spectroscopy [38]. The 
experimental results show a very weak temperature dependence of the exchange split- 
ting in the temperature range 0.3 — 0.85T C , where T c = 1043 K. Our calculations for 
three different temperatures T/T c = 0.3, 0.6 and 0.8 show almost constant d band ex- 
change splitting, in perfect agreement with experiment [38] and previous DMFT(QMC) 
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calculations [7]. 



D. Cr 

The effect of correlations in the case of bcc Cr is manifested through a small en- 
hancement of the density of states at the Fermi level, as shown in Fig. 6. The value of 
electronic specific heat coefficient is 2.88 mJ/K 2 mol, which represents an improvement 
of about 40% relative to the LSDA value of 2.07 mJ/i^mol. Our LSDA+DMFT result 
still underestimates the experimental non-magnetic data of 3.5 mJ/K 2 mol [39] by 15%. 
Although, for the bulk Cr the correlation effects are not very important, one can expect 
strong correlation effects at Cr surfaces, in the light of the observation of essentially 
many-body phenomenon, the orbital Kondo resonance at Cr(001) surface [27]. 

To demonstrate the correlation effect on the real space charge distribution in Fig. 7 
we have plotted the difference of the LSDA+DMFT and LSDA charge densities in the 
bcc (110) plan. As one can see the main effect of the DMFT charge self-consistency is a 
redistribution of charge density, suggesting a suplimentary accumulation of d electrons 
due to correlation effects inside the muffin-tin spheres and a depletion of density in the 
interstitial region. 

E. Fe/Cr multilayer 

In order to test the opportunities of the current LSDA+DMFT approximation fur- 
ther, we have applied this to magnetic heterostructures of alternating magnetic and 
nonmagnetic layers. The most remarkable property of these systems is the giant mag- 
netoresistance (GMR) measured for a parallel/antiparallel configuration of the mag- 
netic moments belonging to the different layer by application of a magnetic field [40]. 
Experimentally was found that the magnetic multilayers grown epitaxially shows an 
enhancement of the electronic contribution to the low temperature specific heat [41]. 
Standard electronic band structure calculations could not reproduce this enhancement 
[42] , which can be an evidence of the correlation effects. Our aim is to check whether the 
correlation effects considered in the EMTO-DMFT approach can lead to an essential 
renormalization of the density of state at the Fermi level N(Ep). 

In the LSDA+DMFT calculations for Fe/Cr multilayer we have chosen different 
values for the average Coulomb interaction and exchange energy for Fe and the two 
types of Cr atoms. These values are listed in Table I. Although this choice is motivated 
by the fact that correlation effects for these atoms could be similar in this structure, 
the present results are more qualitative than quantitative. 

The layer resolved LSDA and LSDA+DMFT densities of states for Fe/Cr multilayer 
are plotted in Fig. 8. Comparing the LSDA and LSDA+DMFT DOS one can see 
that the correlation effects produce a strong peak at the Fermi level. This can give a 
qualitative explanation of the heat capacity data discussed above. 

The calculated LSDA and LSDA+DMFT magnetic moments are shown in Fig. 9. 
Our LSDA result is in agreement with previous electronic structure calculations of Fe 
impurities in a Cr surrounding [43]. We have found that within the LSDA the magnetic 
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multilayer structure consist of ferromagnetically coupled Fe layers with the magnetic 
moments of 1.72/ig per atom. The Cr spacers have very small magnetic moments per 
atom, — 0.05/i£ and 0.09//B, respectively, and they are oriented antiferromagnetically. 
This result is in accordance with previous ab initio study on Fe/Cr superlattices [44]. 
The LSDA Fe magnetic moment is drastically reduced compared with its value in bulk, 
which is attributed to d — d band hybridization between the Fe and Cr states. 

The results of the LSDA+DMFT calculations are essentially different from that in 
LSDA. In particular, the correlation effects result in a slightly increase of the Fe layer 
magnetic moment and a significant polarization of the Cr spacers, see Table I. In the 
case when the correlation effects are considered for all the atoms the selfonsistent cal- 
culation shows the same trend for the magnetic moments: 1.85/ie per Fe atom, 0.42/ie 
per Crl and — 0.19/xs per Cr2 respectively. These results leads us to the conclusion that 
the correlation effects induce a strong polarization on the first Cr layer, which is almost 
non-magnetic according to the LSDA calculations. The appearance of the Cr magnetic 
moments can be attributed to the correlation induced narrowing of the d band, together 
with the Fe-Cr d — d hybridization mechanisms. In the presence of the correlations the 
Fe-Cr d — d hybridization is much stronger than the Cri-Cr2 d — d hybridization. This 
can be seen from majority spins channel of DOS from Fig. 8, where Fe d and Cr x d 
have pronounced peaks at the Fermi level. Due to this significant change in the spin 
up d-chanell of DOS of Fe Fig. 8, the spectral weight is transfered close to Fermi level 
and the effective exchange interaction between the Fe and the first Cr layer is changed 
becoming ferromagnetic. Further investigation of this magnetic ordering as function 
of different parametes U, T is important for the nature of the magnetic coupling in 
multilayer systems. 

Several theoretical approaches have been used to explain the magnetic properties 
of such superlattice structures. Many of these approaches are based on the RKKY- 
like model [45,46], tight-binding models [47] and, recently, on the results of ab initio 
electronic structure calculations [44]. The magnetic coupling studied in the framework 
of these models was shown to result form the interplay between the direct d — d hy- 
bridization of Fe and Cr atoms and indirect exchange through the sp electrons. The 
sp — d coupling [44] was found to be reminiscent of the RKKY interaction only for 
superlattices with more than four Cr layers. 

The calculated electronic specific heat coefficients are listed in Table I. The present 
LSDA+DMFT value of 7.84 mJ/K 2 mol is in good agreement with recent experimental 
study [41]. 

Finally, it is worthwhile to emphasize that the enhanced DOS at the Fermi level, 
having a many-body correlation origin, can play an important role in the GMR, since 
this DOS enhancement is strongly spin-dependent. It is more effective for the majority 
electrons of Fe and Cr! (Oi are the one closer to the Fe layer), giving the result of 
a quasiparticle peak centered at the Fermi level. Our founding is in good agreement 
with the tendency of the enhancement of electronic contribution to the specific heat in 
Fe/Cr magnetic multilayers [41]. 
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V. CONCLUSIONS 



In this paper we present a LSDA+DMFT scheme on the Exact Muffin-Tin Orbitals 
basis set. The present EMTO formalism allows us to combine the many body problem 
with the standard screened KKR(LSDA) method in a self-consistent manner. The 
many-body self-consistency is reached trough self-consistency of self-energy, meanwhile 
the self-consistent charge density is obtained in the conventional LSDA framework. The 
results of our EMTO-DMFT calculation are summarized in Table I, and they are in 
good agreement with former LSDA+DMFT implementations using the LMTO basis 
set. 

Correlations effects in multilayer systems are important for transport properties, 
giving rise to an enhancement of density of states at the Fermi level. We have studied 
a simple magnetic multilayer system, and have shown that the effect of correlation is to 
induce magnetism in non magnetic spacers. We attribute this structure and tempera- 
ture dependent polarization mechanism to simultaneous: (i) Fe-Cr d — d hybridization 
and (ii) narrowing of electronic d bands due to many-body correlation effects. We note 
that the latter mechanism, narrowing of the Cr d band, is significantly stronger than 
the former one. 
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TABLES 





/^LSDA 


/^DMFT 


7LSDA 


7dmft 


fcxpt. 


T 


TJ 


T 

•J 




(Mb) 


(Pb) 


(mJ/K 2 mol) 


(mJ/K 2 mol) 


(mJ/K 2 mol) 


(K) 


(eV) 


(eV) 


Ni 


0.63 


0.42 


5.43 


6.78 


7.02 a 


500 


3 


0.9 


Fe 


2.25 


2.23 


2.43 


2.61 


3.11,3.69 6 


325 


2 


0.9 


Cr 






2.07 


2.88 


3.5 C 


250 


2 


0.9 


Fe(Fe/Cr) 


1.72 


1.75 










2 


0.9 


Cri(Fe/Cr) 


-0.05 


0.15 


6.90* 


7.84* 


8.7±0.7*' rf 


300 





0.0 


Cr 2 (Fe/Cr) 


0.09 


-0.11 













0.0 



a Ref. [33] 
b Ref. [37] 



c Non- magnetic Cr, Ref. [39]. 
d Ref. [41] 

* Values corresponding to the magnetic multilayer. 

TABLE I. Theoretical magnetic moments, /j,, and electronic specific heat coeficients, 7, 

calculated at LSDA and LSDA+DMFT levels. For comparison some experimental electronic 

specific heat coeficients are also listed. The experimental values for 7 include the enhancement 

due to electron-phonon coupling, but this enhancement is not included for the calculated band 

structure values for Fe, Cr and Fe/Cr multilayer. The theoretical 7 values for Ni are corrected 

with 1 + A = 1.24 according to Eq. (19) and Ref. [30]. In the last three columns the parameters 

used in the self-consistent EMTO-DMFT calculations are listed. 
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FIG. 1. The LSDA (dashed line) and LSDA+DMFT (solid line) densities of states for 
fee Ni calculated using the EMTO-DMFT method. A significant reduction of the exchange 
splitting can be evidenced. 
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FIG. 2. Magnetic moment densities of Ni along the nearest-neighbor distance in fee lattice, 
calculated within the LSDA (dashed line) and the LSDA+DMFT (solid line). 
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FIG. 7. The effect of DMFT on the LSDA charge density of bcc Cr. 
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